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Abstract 



We present a novel technique for enhancing Fermi surface (FS) signatures 
in the 2D distribution obtained after the 3D momentum density in a crys- 
tal is projected along a specific direction in momentum space. These results 
are useful for investigating fermiology via high resolution Compton scatter- 
ing and positron annihilation spectroscopies. We focus on the particular case 
of the (110) projection in an fee crystal where the standard approach based 
on the use of the Lock-Crisp- West (LCW) folding theorem fails to give a 
clear FS image due to the strong overlap with FS images obtained through 
projection from higher Brillouin zones. We show how these superposed FS 
images can be disentangled by using a selected set of reciprocal lattice vec- 
tors in the folding process. The applicability of our partial folding scheme 
is illustrated by considering Compton spectra from an Al-3at%Li disordered 
alloy single crystal. For this purpose, high resolution Compton profiles along 
nine directions in the (110) plane were measured. Corresponding highly accu- 
rate theoretical profiles in Al-3at%Li were computed within the local density 
approximation (LDA)-based Korringa-Kohn-Rostoker coherent potential ap- 
proximation (KKR-CPA) first-principles framework. A good level of overall 
accord between theory and experiment is obtained, some expected discrep- 
ancies reflecting electron correlation effects notwithstanding, and the par- 
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tial folding scheme is shown to yield a clear FS image in the (110) plane in 
Al-3%Li. 

PACS numbers: 71.18.+y, 78.70.Ck, 41.60.Ap, 78.70.Bj 
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Typeset using REVT^ 



I. INTRODUCTION 



Recent advances in synchrotron light sources and detector technology have rejuvenated 
interest in high resolution Compton scattering as a probe of fermiology and electron cor- 
relation effects in wide classes of materials. In a Compton scattering experiment one 
measures a directional Compton profile (CP), J{pz), which is related to the twice integrated 
ground-state electron momentum density p{p) by 

JiPz) = jj pip)dp^dpy, (1.1) 

where pz is taken along the scattering vector of the X-rays. In the independent particle 
approximation, 

P(P) = E f AAr)e~'''-'dr\ (1.2) 
b,k '' 

where ipt.ki'^) is the electron wavefunction in band b and state k. The summation in Eq. 



( |1.2| ) extends over all initial states that are excited in the scattering process. The CP thus 
contains signatures of the Fermi surface (FS) breaks and correlation effects in the underlying 
three-dimensional (3D) momentum density p(p). 

The determination of FS geometry via the Compton spectroscopy is complicated by the 
presence of the double integral in Eq. ( |1.1[ ). As a result, FS breaks in p{p) do not in general 
induce rapid variations in J{pz), but only in the first derivative of the profile. Moreover, 
breaks in p(p) caused by FS crossings are scattered throughout momentum space with 
weights given by the appropriate matrix element. Therefore, if this matrix element happens 
to be small in a specific case,lli then the related FS feature in the CP will be intrinsically 
weak and difficult to measure. The standard approach to tackle these problems is to measure 
CP's along many different directions, and invert Eq. ( p. . 1|) to obtain p{p) by using one of 
the available reconstruction methods.lll'lll The matrix element effects in p{p) can then be 
removed in principle by constructing the occupation number density in p space, N{p), via 
the so-called Lock-Crisp- West (LCW) folding procedureEl; 
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N{p) ^Y.PiP + G) = T.i: ^ (1-3) 

G G b,k 

where the sum over G runs over all reciprocal lattice vectors. The right hand side of Eq. 
( p..3[ ) represents the occupation number density in the repeated zone scheme.lii 

In practical applications, an accurate reconstruction of p{p) along the preceding lines re- 
quires the demanding task of making CP measurements along a suitable group of directions. 
However, in many cases the 2D projection of the momentum density, p{px,Py) = J p{p)dpz, 
and the corresponding 2D occupation number density, N{px,Py), along a judiciously chosen 
direction can give most of the important features of the FS. As shown in Sec. || it is then 
sufficient to measure a set of CP's with scattering vectors all lying in the plane of projection, 
yielding a significant reduction of measuring time. The (110) plane in cubic lattices encom- 
passes many of the high symmetry points and lines in the Brillouin zone and, therefore, a 
cross-section in this plane captures much of the character of the underlying 3D Fermi surface. 
In this vein, the (110) projection has been exploited in a number of Compton scattering and 
positron annihilation investigations&^^ on bcc crystals where the (110) projection turns out 
to be reasonably free of complications of overlapping FS images from higher Brillouin zones. 
In sharp contrast, it has not been possible to use this approach in fee crystals due to overlap 
problems which often distort the FS image severely. 

With this motivation, we have examined the question of how the aforementioned prob- 
lems of overlap may be avoided in a projected 2D momentum distribution. Our analysis 
reveals that the problem can be solved by using a selected subset of reciprocal vectors in the 
folding procedure of Eq. (pT^ ). To our knowledge, such an approach has not been presented 
previously in the literature. Although we focus in this article explicitly on the specific case 
of the (110) projection in an fee crystal, these ideas are generalized straightforwardly to 
other projections in cubic lattices. The treatment of projections in general lattices will be 
taken up elsewhere. 

We illustrate the applicability of our novel folding procedure by considering Compton 
spectra from an Al-3%Li disordered alloy (fee) single crystal specimen. For this purpose, 
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nine high resolution CP's have been measured along a set of directions in the (110) plane, and 
corresponding highly accurate theoretical profiles have been computed within the KKR-CPA 
first-principles framework. Intercomparisons between theory and experiment, and between 
results based on the use of a selected and a full set of reciprocal lattice vectors in the folding 
process, clearly show the value of our partial folding approach in enhancing FS signatures 
in the data. 

Some comment on our choice of the Al-3%Li disordered alloy for this study is appropriate. 
Al, which possesses an almost spherical FS (in the extended zone), is an interesting base 
metal for exploring the evolution of fermiology with changes in the electrons/atom ratio via 
alloying on the one hand, and for investigating electron correlation effects on the momentum 
density in multivalent systems on the other. Metallurgically, Al is quite "exclusive" in that 
only a few percent impurities can be dissolved in the solid solution a phase.@ For example, 
in Al-Li the a phase extends at room temperature to only about 3 at% Li concentration. 
Al-Li alloys, with the possible addition of other elements, are also of great technological 
interest due to their high strength-to-weight ratio.B For these fundamental and practical 
reasons, the electronic structure and the nature of bonding in Al-Li alloys have drawn 
considerable interest .il^ Insofar as the Compton technique is concerned, relatively little 
by the way of high resolution work on disordered alloys is available in the literature despite 
the fact that, unlike traditional spectroscopies which rely on quantum oscillatory phenomena, 
Compton scattering is not affected by short electron mean free paths and is thus particularly 
suited for investigating disordered alloys. For reasons outlined in this paragraph, there is a 
strong motivation for undertaking a state-of-the-art Compton study of the Al-Li system as 
a means of gaining insight into the electronic structure of the system and the possible role 
of fermiology in destabilizing the a phase of Al.S 

This article is organized as follows. Section y briefly delineates some reconstruction 
aspects needed later. It is shown that it suffices to measure a set of CP's with scattering 
vectors all lying in one plane to obtain the 2D projection of the momentum density onto 
that plane. In Sec. |T|, the new folding procedure involving selected reciprocal lattice vec- 

6 



tors is formulated. Section |I^ gives pertinent details of the Compton experiments, while 
Sec. |V| briefly describes the relevant computational details involved in obtaining theoretical 
momentum densities and Compton profiles. Section |VT| presents and discusses the new the- 
oretical and experimental Compton results as well as the applicability of our partial folding 
procedure to these data. Some concluding remarks are collected in Sec. |V14 

II. RECONSTRUCTION 

The 3D momentum density can be reconstructed by using the Fourier-Bessel method in 
which a so-called B{r) function is defined by 

B{r) = J p{p) exp{—ip ■ r)dp. (2.1) 

It follows from Eqs. (^) and (^) that 

B{0,0,z) = J J{p,)exp{-tp,z)dp,. (2.2) 

Thus, from a number of measured CP's, corresponding i?(r)'s are obtained on radii in r 
space given by the directions of the respective scattering vectors. Next, a fine mesh is set up 
in r space and B{r) is obtained at every mesh point by interpolation. This makes it then 
possible to carry out the inverse Fourier transform 

p{p) = {2n)~^ j B{r) exp(ip ■ r)dr (2.3) 



to obtain p{p). We refer the reader to Ref. |T7|for a detailed discussion of our reconstruction 
scheme. 

Here we briefly show how a 2D projection of the momentum density along a given 
direction can be reconstructed from Compton measurements in which all the scattering 
vectors lie in the perpendicular plane. Without loss of generality, we take [001] to be the 
projection direction. From the inverse Fourier transform Eq. ( |2.3| ) the projection onto the 
(001) plane is: 
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p{Px,Py) = J p{p)dpz = (27r) ^ J dpz J J J B{x,y,z) 

X exp[i{pxX + PyU + Pzz)]dxdydz. (2.4) 

The integration with respect to Pz produces 2tt6{z), yielding 

= (27r)"2y" J B{x,y,0)exp[i{p^x+pyy)]dxdy. (2.5) 

But, as seen from Eq. ( p.2|) above, B{x,y,0) can be obtained by interpolation on a fine 
mesh in the (x, y) plane from measurements in which all scattering vectors lie in that plane. 

III. FORMULATION OF PARTIAL FOLDING TECHNIQUE 

This section delineates briefiy how the (110) projection of the LCW-folded momentum 
density in Eq. ( [1.3|) can be divided for an fee crystal into two identical patterns where one 
pattern is shifted by 27r/a along the [001] direction with respect to the other. 

We consider an fee lattice with primitive translation vectors: 

a=^{x + y) ; b=^{y + z) ; c = ^{z + x) (3.1) 
where x, y,z are unit vectors along the Cartesian axes. A general lattice vector is 

R = ua + vb + wc ; u,v,w are integers. (3.2) 

Similarly, the primitive translation vectors of the (bcc) reciprocal lattice are : 

27r ,^ ^ ^, „ 27r , ^ ^ ^, 

A = —{x + y- z) ; B = —{-x + y + z) ; 
a a 

C = —(x-y + z), 

a 

and the general reciprocal lattice vector is 

G = hA + kB + IC ; h,kj are integers. (3.3) 

Because we are interested specifically in projecting onto the (110) plane, it is convenient to 
work with new unit vectors s and t which lie along and perpendicular to [110]: 
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and thus 



with 



while 



with 



R = RsS + Rtt + R,z (3.5) 



Rs = ^^i2u + V + w) ; Rt = ^-^{-v + w) ; 

Rz = ^{v + w), (3.6) 



G = GsS + Gtt + G,z (3.7) 



Gs = —V2h ; Gt = —V2{-k + l) ; 
a a 

2-71 

G, = —{-h + k + l). (3.8) 
a 



Likewise, we have 



p = PsS + pft + and r = rgS + rtt + r^z. 
The momentum density projected onto the (110) plane is 

J p{p)dPs = p{pt,Pz)- (3.9) 
Having established the relevant notation, we write the momentum density Eq. ( |1.2|) as 

p{p) = p{Ps,PuP.) = Y.II drdr' ^.,fc(r)^,%(r')e-''(^-^'). (3.10) 

, iJ J crystal 



b,k 



I crystal 

Writing the wavefunctions in the Bloch form and converting the integral over r in Eq. ( [L.2| 
into a sum of integrals over a single cell one obtains 



p{Pt,Pz) = dpsY^Y. / / drdr'ub^k{r)ul,,{r') 

b,k R,R' 

^gifc-(H+r)g--ifc-(R'+r')g-jp(fi-fi')g-*p(r--r')_ (3-11) 



The two-dimensional projection of the occupation number density N{pt,Pz) obtained by 
LCW-folding over all projected reciprocal lattice vectors is 

Nipt,P.) = E PiPt + Gt,P. + G,) (3.12) 

Gt,Gz 

where the sum runs over all sets {Gt, G^)- 

The formulation so far leading to Eq. ( p.l2|) has been quite general. The nature of the 



problem as well as its solution may be delineated with reference to Fig. |l] which shows 
schematically the projection of the bcc reciprocal lattice onto the (110) plane. The lattice 
points in the (110) plane going through the origin (for which h = in Eq. (3.8)) lie at 
the positions of the circles. The two nearest (110) planes {h = ±1) lie at a distance of 
Gs = {2n/a)\/2 in front of and behind this central plane and the corresponding lattice 
points project into the crosses. The next pair of planes projects to the circles and then the 
process repeats itself. The entire bcc reciprocal lattice is thus divided into two sublattices 
related to the circles and the crosses with associated 2D Brillouin zones shown by light and 
dark lines respectively.0 It is easily seen now that when the images of the Fermi surface 
given by the periodic LCW-folded 3D-density are projected along [110], the contributions 
from the circled and crossed sublattices will overlap and get entangled. However, since the 
circled and crossed sublattices are related by a simple shift of 271 /a along the z axis, if 
the contribution of one of these sublattices can be suppressed in the folding procedure, the 
projection will no longer suffer from the aforementioned complication. 

The preceding strategy may be carried out more formally as follows. In the (s, t, z) space, 
Eqs. ( p.8|) give the coordinates of various bcc lattice points generated via the integers /i, k 
and /. Different layers perpendicular to the [110] direction correspond to different values of 
h. It is easily shown then that the circled sublattice in Fig. |I] involves only even values of 
the index h where Gt and Gz (ignoring prefactors in Eqs. ( |3.8| )) are either both even or both 
odd. Similarly, the crossed sublattice is given by odd values of h where Gt and Gz (again, 
ignoring prefactors) possess opposite parities. 

Bearing these considerations in mind, we divide the summation in Eq. ( p.l2|) into two 
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sums involving even and odd values of h and we use Eq. This yields 



hieven h:odd 

Gt,Gz Gt,Gz 

h:even „ r r 

Gt,G^ b,k R,R' •' 

-ips{Rs-R's+rs-r',) -i{pt+Gt){Rt-R't+rt-r'^) -i{p,+Gz){R:,-R'z+r,-r'^) 



+ a similar expression with h odd. (3.13) 
We now make use of the following identities: 

I c?p,e-^f=(^-^'=+'--'-^) = Sn'^nJir, - r^) (3.14) 

J2 ^-iGt{n-r',)-iG4r.~K) = ^(^^ _ _ (3,15) 

Gt,Gz 

Y^Y,^^(M-p^)R^^^(k.~p.)R. ^ ^^^^^^^ + G^)6{K,p, + G^) , (3.16) 

Rt Rz 

where for the sake of simplicity we ignore forefactors, 5(a, h) denotes the Kronecker delta 5afi, 
and (G", G^) are the components of an arbitrary reciprocal lattice vector G°'.iil Furthermore, 

■r. ,r, ^,r^ / n ™^ I 1 foreVCU/l 

^-iGt{Rt-R't)-tGz{Rz-R'z) ^ J (3-17) 

^-^1{RZ-K) fQJ. 



The last identity is easily seen by using Eqs. ( |3.6|) and ( pTSf) and realizing that v + w = 
{2/a)Rz- The final result is 

N'''''''{pt,P.)=c Y.^{h,pt + Gl)6{Kp. + Gl) (3.18) 

G?,G| b,k 

and 

h-.odd c\ 

N'"'\pup^) = c J2 Y^ikt,pt + G",)5{k,,p, + — + G:), (3.19) 

G°,G° b,k ^ 

where the constant c collects the irrelevant prefactors and Gf^°^ and Gf-°^ are the components 
of G with even (odd) 

The right hand side of Eq. ( ^.18] ) is just the projected two-dimensional k space occupation 
sampling function in the repeated zone scheme for the set of reciprocal lattice vectors with 
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even h in Eq. ( |3.3| ), while the right hand side of Eq. (|3.19|) fulfills the same function for the 
set with odd h. Thus, the set of Fermi surface images appearing in Eq. (|3.18| ), reappear in 
Eq. ( |3.19| ) but shifted by 27r/a along z. By carrying out the LOW folding over only one of 
these two sets of reciprocal lattice vectors, the Fermi surfaces corresponding to the other set 
are removed and thus a picture is obtained which is not complicated by superposition. 

Since N^^^"^ and A^"'^'^ are identical (apart from the shift) and their sum yields the occu- 
pation function N{pt,Pz), it follows that N^^'^^ and A^"'^'^ each carry half the intensity of the 
occupation function. 

Finally, it is clear that the present procedure can be used straightforwardly in the analysis 
of positron annihilation 2D-ACAR spectra from an fee crystal. 



IV. EXPERIMENT 

A single crystal of Al-3at%Li was grown from the melt by the Bridgman-Stockbarger 
method under a pressurized argon atmosphere to reduce loss of Li. The composition was 
determined by both Atomic Absorption Spectrometry and Proton Induced Gamma Emission 
and the final specimen, 10 mm in diameter and 1.5 mm thick, was found to contain 3.0 at% 
Li. The Compton profiles along nine directions, all lying in the (110) plane (see Fig. 
were measured with the KEK-AR spectrometer.^! The incident photon energy was 60 keV 
which is high enough so that the impulse approximation would be satisfied quite well. The 
overall momentum resolution in the measurements is estimated to be 0.12 a.u. The total 
number of accumulated counts under the profiles was about 1 X 10^ The necessary energy 
dependent corrections for absorption, detector and analyzer efficiency, and scattering cross- 
section were applied. The contribution of double scattering was simulated via the Monte 
Carlo program of Sakai.0 The integrated intensity of the double scattering events was found 
to be 10 % of that of the single scattering events. By using Eq. ( p.2| ), the nine measured 
profiles yield the values of B{z) along nine radial directions all of which lie in the (110) 
plane in real space. B{r) was then interpolated over an area of 50 x 50 a.u.^ on a fine square 
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mesh of 500 x 500 r points in the (110) plane. The projection p{pt,Pz) of the 3D momentum 
density was obtained via an inverse 2D- Fourier transformation of the interpolated B{r). 
The p{pt,Pz) thus obtained was partially folded following the procedure of Sec. ^ to obtain 
the 2D-occupation number density. 



V. COMPUTATIONS 

The computation of the electronic structure of Al-3at%Li was carried out within the all 
electron charge self-consistent KKR-CPA framework and is parameter-free. A lattice con- 
stant of a = 7.6534 a.u. was assumed. The underlying KKR-CPA methodology is described 
in Refs. |35|-p9|. The relevant Green's function formulation for treating the momentum density 



and Compton profile in disordered alloys is given in Refs. ^ and ^ Exchange and correla- 
tion effects were incorporated within the von Barth-Hedin local density approximation.Eijli 
In order to obtain the Compton profile the momentum density p{p) was first evaluated over 
a fine mesh of 48 x 4851 x 1777 p points covering momenta up to Pmax ~ 5 a.u. Compton 
profiles with scattering vectors along the same directions as the experimental ones were then 
computed accurately by evaluating the integral of Eq. ( |1 . 1| ) . The accuracy of the computed 
profiles is about 1 part in 10^. The nine theoretical profiles were convoluted with the exper- 
imental resolution and then treated in a manner identical to that of the nine experimental 
profiles. 



VI. APPLICATION TO AL-3%LI ALLOY 

We briefiy compare first the theoretical and experimental directional profiles by taking 
the [100] CP as an example; results along other directions are similar. After subtraction 
of the core contribution the two profiles of Fig. ^ are area-normalized to yield the same 
number of electrons over the range ±4 a.u. The theoretical CP is seen to be higher than 
the experimental one at low momenta with the situation reversing itself at high momenta. 
The Fermi cut-off in the first derivative of the CP around 1 a.u. is sharper in the theory; 
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this effect is also evident in the second derivative where the corresponding calculated peak 
is higher and narrower compared to the measurements even though the theoretical CP's 
and their derivatives include resolution broadening. Similar discrepancies between theory 
and experiment have been observed previously in other systems and can be ascribed to 
the failure of the present LDA-based independent particle model to properly account for 
correlations in the electron gas.iHii These correlation effects result in the shift of spectral 
weight from below to above the Fermi momentum (pp) and a reduction in the size of the 
break (Zk) in the momentum density at pp- Concerning fine structure, features around 0.3 
a.u. and 0.6 a.u. in the first derivative of the experimental CP are reproduced reasonably 
by theory. These features are quite similar to those found in Al in the positron annihilation 
(ID-angular correlation) experiments^ as well as in the recent high-resolution Compton 
scattering measurements.^ In the case of Al, these have been explained by the fact that the 
Fermi surface of Al deviates from a sphere (in the extended zone) and bulges on the hexagonal 
zone face near the W-K-W and W-U-W zone edges;S this mechanism is presumably at play 
in Al-3%Li as well since the addition of 3% Li induces relatively little change in the Fermi 
surface of Al. 

Before discussing the effects of folding, it is helpful to consider Fig. ^ which shows 
contour maps of the reconstructed experimental and theoretical momentum density in the 
(110) plane. We emphasize that an integration along [110] (i.e., Ps) is implicit here as well as 
throughout this article. Therefore, what is shown in all cases is the (110) 2D projection of 
the underlying 3D momentum density distribution. As expected, the data of Fig. ^ display 
two- fold crystalline symmetry, but it is difficult to ascertain the details of the FS, except 
perhaps that it is roughly spherical in shape (in the extended zone) with some bulging along 
the (110) directions. 

Insight into the nature of our partial folding scheme is provided by Fig. |^ which considers 
the case of the free electron model where the Fermi sphere corresponds to 2.94 electrons/atom 
appropriate to the number of conduction electrons in Al-3%Li. Figure |^(a) has been obtained 
via the conventional LCW-folding procedure in which all reciprocal lattice vectors are used, 
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whereas Fig. ^(b) includes the contribution of only one of the sublattices as discussed in 
Sec. |T|. The imprint of the free electron FS (repeated periodically) is obvious in the 2D- 
occupation number density of Fig. ^(b), even though this is far from being the case in Fig. 
^(a) which shows the appearance of four cross-like features around the L symmetry points 
due to overlapping Fermi surfaces. Furthermore, complicated structures arise in Fig. |^(a) 
around the zone center. It is striking how the two overlapping Fermi surfaces in Fig. |^(a) 
are disentangled in Fig. ^(b).0 These results show clearly the power of the present partial 
folding approach for analyzing FS images in terms of the (llO)-projection in an fee crystal. 

The application to Al-3%Li is considered in Fig. ^ which compares theoretical and ex- 
perimental (llO)-projected, partially folded, 2D occupation densities. On the whole, the FS 
contours in Fig. ^ are rather like those in Fig. ^(b), indicating that the FS of the alloy is 
more or less spherical, despite some deviations. Note that the FS contours in Fig. |^ are not 
symmetrical with respect to the U-K lines, even though this is strictly true for the free elec- 
tron case in Fig. ^(b). This is a consequence of the asphericity of the FS; the FS bulges more 
outside the hexagonal zone face near the W-K-W zone edge (where two hexagonal faces of 
the first zone intersect) than near the W-U-W zone edge (intersection between a hexagonal 
and a square face of the first zone). This can also be described as a slight flattening of the 
FS in the vicinity of U, combined with a small expansion near K. Theory and experiment 
are in good overall accord in Fig. although the fact that the oblong maximum around 
the L point in Fig. ^(b) is slightly more rotated compared to Fig. |^(a) suggests that the 
aforementioned asphericity of the FS is somewhat stronger in the experiment. 

Finally, we consider radial cuts through the 2D momentum distribution of Fig. ^; a typical 
cut (along F-L) is shown in Fig. |^. It is well known that the value of the Fermi momentum 
Pf is not given correctly by the inflection point in the profile of Fig. 0(a), or equivalently, 
the position of the minimum in the first derivative of Fig. |^(b).@ The situation may be 
simulated by considering the free electron case. Here the 2D-distribution corresponding to 
the profile of Fig. |^(a) is a semi-circle and the first derivative is negative and diverges at pp. 
When the semi-circle is convoluted with the experimental resolution, a tail appears beyond 
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Pp and the divergence of the first derivative &t pp disappears. The inflection point defining 
the minimum of the first derivative then moves below pp to an extent which depends on 
the width of the resolution function. Accordingly, by using a free electron semi-circle of 
radius 0.919 a.u. and the present experimental resolution of 0.12 a.u., we have determined 
that the position of the inflection points in cuts such as those of Fig. |^ must be increased 
by 4.2 % in order to obtain the correct value of pp. Applying this correction, theoretical 
and experimental pp values in the (110) plane can be mapped out from the data of Fig. 

The results of the analysis are given in Fig. and show an excellent level of accord 
between theory and experiment. Nevertheless, a careful examination of Fig. reveals a 
slight contraction of the experimental FS (compared to theory) near U in combination with 
a slight expansion near K, consistent with our earlier discussion in connection with Fig. ^ 
above. 



VII. SUMMARY AND CONCLUSIONS 

We address the question of how Fermi surface (FS) signatures can be enhanced in the 2D 
distribution obtained by projecting the 3D momentum density along the [110] direction in 
an fee crystal. The standard LCW folding procedure invoked in this connection involves a 
summation of the momentum density over all reciprocal lattice vectors to obtain the electron 
occupation number density in the system. This procedure has been used to produce a rea- 
sonably clear FS imprint, even when the (110) projected momentum density is considered, 
in several studies of bcc crystals, but it has not been successful in fee crystals. We show, 
however, that the (110) projection of the LCW-folded momentum density in the fee lattice 
can be viewed as a superposition of two equivalent FS images which are shifted by 27r/a 
along the [001] axis, and that these two images can be disentangled if one uses a selected 
subset of reciprocal lattice vectors in the folding process. Our novel folding procedure is 
particularly well suited for investigating fermiology-related issues via high resolution Comp- 
ton experiments where a full 3D reconstruction of the data is far more demanding than a 
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2D-reconstruction using measured directional profiles all of which involve scattering vectors 
in the (110) plane. The technique will also be valuable in analyzing positron annihilation 
(2D-ACAR) experiments in fee lattices. 

An extensive application of the aforementioned partial folding procedure to the case of 
an fee disordered alloy is presented. To this end, we have measured nine high resolution di- 
rectional Compton profiles on a single crystal specimen of Al-3at%Li with scattering vectors 
in the (110) plane. We have also carried out highly accurate computations of corresponding 
theoretical profiles within the self-consistent, parameter free KKR-CPA framework. Occu- 
pation number densities based on theoretical as well as the experimental profiles obtained 
via partial folding display clear images of the FS. A reasonable level of accord is found be- 
tween theory and experiment with respect to the overall shape of the directional profiles as 
well as the fine structure in the first and the second derivatives. The experimental results 
with regard to the FS are also in accord with KKR-CPA predictions, although there is an 
indication of a slightly stronger bulging of the measured Fermi surface outside the hexagonal 
zone face along the W-K-W than along the W-U-W zone edge. Upon adding 3% Li, the 
FS of Al remains essentially free electron hke, some overall shrinking of its size due to the 
reduced e/a ratio notwithstanding. Large FS sheets in Al thus show no sign of an anomaly 
which might destabilize the a phase. A Fermi surface driven mechanism in this connection 
will therefore need to focus on small features in the FS of Al (electron pockets in the third 
conduction band) which can undergo substantial change and may even disappear with de- 
creasing electron concentration or altered electron lattice interaction upon alloying. 
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FIGURES 

FIG. 1. The bcc reciprocal lattice points (corresponding to an fee crystal) are shown projected 
onto the (110) plane. 2 denotes the [001] and t the [110] direction. As discussed in the text, open 
circles refer to points with even h in Eq. (3.8) and crosses to points with odd h. Outlines of the 
Brillouin zones for the sublattices of circles and crosses are shown by thin and thick solid lines 
respectively. 

FIG. 2. The nine directions in the (110) plane along which the Compton profiles were measured 
and computed in this study. 

FIG. 3. Experimental and theoretical directional Compton profiles in Al-3%Li and the asso- 
ciated first and second derivatives are compared along the [001] direction. The theoretical results 
have been broadened to reflect the experimental resolution. 

FIG. 4. Contour maps of the reconstructed experimental and theoretical momentum densities 
in the (110) or {pz,Pt) plane. Note that the results shown represent the 2D-distribution obtained 
after the 3D momentum density is integrated along the [110] direction. The first Brillouin zone 
boundary and the high symmetry points in the (110) plane are indicated. 

FIG. 5. Contour maps of the 2D occupation number density (after projection along [110]), 
N{kz,kt), in the (110) plane for the free-electron model at the electron concentration of Al-3%Li. 
(a) Conventional LCW folding using all reciprocal lattice vectors, (b) Present partial LCW folding 
scheme using selected reciprocal lattice vectors. Experimental resolution broadening is included in 
the calculations. The Brillouin zone boundary and the projected free-electron sphere are shown. 
+ signs indicate high density regions. 

FIG. 6. Contour maps of the theoretical and experimental 2D occupation number densities 
(after projection along [110]), N{kz, kt), in the (110) plane based on the present partial LCW-folding 
scheme in Al-3%Li. Resolution broadening is included in the theory. See caption to Fig. |5| for 
other notational details. 
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FIG. 7. Radial sections along the F-L direction through the experimental and theoretical 
[110]-projected momentum distributions of Fig. ^ are shown together with the corresponding first 
derivatives. 

FIG. 8. Fermi momentum pp in the (110) plane in Al-3%Li disordered alloy. Open circles 
give results deduced from the present Compton profile measurements, while the KKR-CPA based 
theoretical predictions are given by filled circles. The solid circular arc indicates the free electron 
momentum of 0.919 a.u. for Al-3%Li. 
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